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Abstract 



The generating functionals (GF) method in Bogolyubov's formulation and its 
application for particle physics is considered. Effectiveness of the method is illus- 
trated by two examples. So, GF method can be used as the technical trick solving 
the infinite sequence of algebraic equations. We will consider the example, where 
• GF allows express the multiplicity distributions (topological cross sections) through 

the particles correlation functions (inclusive cross sections) to 'predict' so called the 
Koba-Nielsen-Olesen scaling. We will use the GF to define validity of the thermal 
description of the multiple production phenomena also. It will be seen that this will 
lead to the 'correlations relaxation condition' of N.N.Bogolyubov. This will allow 
Qh| to offer the experimentally measurable criteria of applicability of thermodynamical 

Q_i' description of multiple production processes. In results we will find the closed form 

of perturbation theory applicable for kinetic phase of nonequilibrium processes. It 
is shown a way as the approach may be adapted to the definite external conditions. 

> ■ 
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1 Introduction 

It is hard to imagine modern particles physics without such fundamental notions as, 
for instance, the phase transitions, topological defects, taken from statistical physics. 
This extremely fruitful connection among two branches of physics based on the euclidean 
postulate [1]: the formulae of particle physics are coincide with corresponding formulae 
of statistical physics if the transformation t — ^ it is applied. But this coincidence exist 
iff the media is equilibrium only, since the time order of physical process becomes lost 
after the transition to imaginary time it. So, the particles static properties only can be 
considered by euclidean field theories. 

The euclidean postulate does not 'work' for arbitrary element of S'-matrix and, by this 
reason, there is not, at first glance, general connection between particles and statistical 
physics. Our aim is demonstrate this connection considering the multiple production 
example, staying in the real-time theory frame. 

The multiple production is a typical dissipative process of the incident kinetic energies 
transition into the energies (masses) of produced particles. This is the nonequilibrium 
process and the fluctuations, generally speaking, may be high in it. Experimental data 
confirms this general expectation at the mean multiplicities region, when n ~ n [2] . 
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Considering multiple production we would like to note firstly that the mean multiplic- 
ity n of hadrons for modern accelerator energies (~10 Tev) is large n(s) ~100. So, it is 
practically impossible to describe the system with iV = 3n — 10 ~ 300 degrees of freedom 
using ordinary methods. 

Secondly, it is natural to assume that the entropy S tends to maximum with rising 
multiplicity n and reach the maximum at n ~ n max ~ y/s, since the dissipation take place 
in the vacuum (presumably with zero energy density) 2 . But the experiment shows that 
at high energies n ~ n(s) ~ In 2 s « n max (s) are essential. This means that there is not 
total dissipation of incident energy in considered thermalization process [3]. Absence of 
thermalization may be a consequence of hidden conservation laws [4] . 

We would like to adopt following fundamental principle of nonequilibrium statistics 
introduced by N.N.Bogolyubov [5]. It is natural to assume that the system evaluate to 
the equilibrium in such a way that the 'nonequilibrium' fluctuations in it should tend to 
zero. In the frame of Bogolyubov's principle the quantitative measure of 'nonequilibrium' 
fluctuation is the mean value of correlation functions and, therefore, this quantities should 
tend to zero when the media tends to equilibrium. 

In our interpretation the Bogolyubov's correlations relaxation principle means follow- 
ing. So, for nonequilibrium state presence 'nonequilibrium' fluctuations in the form of the 
macroscopic flow of, for instance, energy e is natural. Then the mean value of m-point 
correlation functions K rn can not be small as the consequence of macroscopic flow. But 
in vicinity of equilibrium the macroscopic flows should relax and, accordingly, the mean 
value of correlation functions should be small, K m 0. To characterize the equilibrium 
one may consider also the particles, charge, spin, etc. densities macroscopic flows and 
theirs relaxation. 

We would like to show in result that the correlations relaxation principle leads to the 
quantitative connection with real time thermodynamics of Schwinger-Keldysh type 3 [6]. 
Just for this purpose the generating functionals (GF) method of Bogolyubov will be used 
since it allows to find the quantitative connections, where the euclidean postulate does 
not applicable. 

We will use more natural for particles physics microcanonical formalism. In this for- 
malism the thermodynamical 'rough' variables are introduced as the Lagrange multipliers 
of corresponding conservation laws. Theirs physical meaning are defined by correspond- 
ing equations of state. So, if the fluctuations in vicinity of solutions of corresponding 
equations are Gaussian then one can use this variables for description of the system. 
Corresponding condition is the Bogolyubov's correlations relaxation condition. 

Formally, the generating functions method presents the integral transformation to new 
variables. One can choose them as the 'rough' thermodynamical variables. To describe 
the far from equilibrium system we will introduce the 'local equilibrium hypothesis'. In 
its frame the preequilibrium state consist from equilibrium domains. In this case new 
variables should depend on the coordinates of domain and, in result, we are forced to use 

2 This consideration lie in the basis of earliest Fermi-Landau 'statistical' model of hadrons multiple 
production. 

3 Last one includes the nonequilibrium thermodynamics also. 



the generating functional (GF) formalism. 

We will consider two example to illustrate effectiveness of the GF method. In Sec. 2 
we will consider the transformation (multiplicity n — > activity z) to show the origin of 
the Koba-Nielsen-Olesen scaling (KNO-scaling) 4 . 

In Sec. 3 we will investigate a possibility of temperature description of the multiple 
production processes. We will consider for this purpose the transformation (particles 
energies e — > temperature 1//3) to find the S'-matrix interpretation of thermodynamics. 
It will be shown that this interpretation would be rightful iff the correlations are relax. 

In Sec. 4 we will use this interpretation to formulate the perturbation theory in the 
case when (5 and z are local coordinates of temperature (x,t) [7]. One can use this closed 
form of perturbation theory for description of nonequilibrium media (in kinetic phase) 
and for description of the multiple production process as well. 



2 KNO-scaling 

We would like start from note that the generating functions method allows connect in- 
clusive spectra fa [8] and exclusive cross sections & n (s). One can use for this purpose the 
normalization condition: 

n\ _ 

n=k {n-k)\ 

where, as usual, 

du k (q) = f[d%/(27r) 3 2e( qi ), e(q) = ^/q 2 + m 2 
i=i 

, is the Lorentz-covariant element of phase space. 

Eq.(2.1) can be considered as the set of coupled equations for a n . One may multiply 
both sides of (2.1) on (z — l) k /k\ and sum over k to solve them. We will see that this is 
equivalent of introduction of 'big partition function' S(z), where z is the 'activity': the 
chemical potential ~ \nz. 

We will find in result of summation over k that 

k K - n 

Then, assuming that S(z) is known, 




where the closed contour C includes point z — 0. Here S(z) is defined by left hand side 
of (2.2) and is the generating function of a n . 

4 In privet discussion with one of the authors (A.S.) at summer of 1973 Z.Koba noted that the main 
reason of investigation leading to the KNO-scaling was just the GF method of N.N.Bogolyubov. 



(2.2) 



The coefficients C m in decomposition: 



m-^HE^-^O, (2.4) 

are the (binomial) correlators. Indeed, 

d = h =n, C 2 = f 2 - {A} 2 , C 3 = / 3 - 3/ 2 {/i} 2 + 2 {A} 3 (2.5) 
an so on. If C m = 0, m > 1, then <r n is described by Poisson formulae: 

o n = a tot e- n{ -^-. (2.6) 
nl 

It corresponds to the case of absence of correlations. 
Let us consider more week assumption: 

C m (s) =7m(Ci(s)) ro , (2.7) 

where 7 m is the energy independent constant. Then 

hiE(z,8)=?:%{(z-i)n(8)r. (2.8) 

m=l m - 

To find consequences of this assumption let us find the mostly probable values of z. The 
equation: 

d 

n = z—lnZ(z,s) (2.9) 

oz 

has increasing with n solutions z(n, s) since E(z, s) is the increasing function of z, iff 
S(z, s) is the nonsingular at finite z function. Last condition has deep physical meaning 
and practically assumes that absence of first order phase transition [9]. 
Let us introduce new variable: 

\ = (z-l)n(s). (2.10) 
Corresponding eq.(2.9) looks as follows: 

" 'l + ^Ul^A). (2.11) 



n(s) \ n(s) J d\ 

So, with 0(X/n(s)) accuracy, one can assume that 

A ~ X c (n/n(s)). (2.12) 

are essential. It follows from this estimation that such scaling dependence is rightful at 
least in the neighborhood of z — 1, i.e. in vicinity of main contributions into o to t- This 
gives: 

n(s)a n (s) = <j tot (s)ip(n/n(s)), (2.13) 



where 

i){n/n{s)) ~ E(X c (n/n(s))) exp{n/n(s)\ c (n/n(s))} < 0(e~ n ) (2.14) 

is the unknown function. The asymptotic estimation follows from the fact that A c = 
X c (n/n(s)) should be, as follows from nonsingularity of S(z), nondecreasing function of 
n. 

The estimation (2.12) is rightful at least at s — > oo. The range validity of n, where 
solution of (2.12) is acceptable depends from exact form of S(z). Indeed, if lnS(z) ~ 
exp{7A(z)}, 7 = const > 0, then (2.12) is rightful at all values of n and it is enough to 
have the condition s — > oo. But if InE(z) ~ (1 + aA(z)) 7 , 7 = const > 0, then (2.12) is 
acceptable iff n « n 2 (s). 

Representation (2.13) shows that just n(s) is the natural scale of multiplicity n [10]. 
This representation was offered firstly as a reaction on the so called Feynman scaling for 
inclusive cross section: 

k x 

fk(qi, Q2, Qk) ~ II —, r- (2.15) 

i=l e \1i) 

As follows from estimation (2.14), the limiting KNO prediction assumes that a n = 
0(e~ n ). In this regime S(z, s) should be singular at z — z c (s) > 1. The normalization 
condition 

dZ(z,s) _ 
-^-U = n(«) 

gives: -2 C ( S ) = 1 + j/n(s), where 7 > is the constant. Note, such behavior of big par- 
tition function E(z, s) is natural for stationary Markovian processes described by logistic 
equations [11]. In the field theory such equation describes the QCD jets [12]. 

It is known that at Tevatron energies the mean hadrons multiplicity rise with trans- 
verse momentum. The associated mean multiplicity is 

Ci{q tr ) = n(q tr ) - 



J2 n dcr n /dq tr ' 
So, if 



k 



C m (q tr ) = 7m(Ci(%r)) m : /fc(gi,92,-,5fc) ~ H -f^:^(Qtr) , 

i=l £ {Qi) 

then: 

_ dajdqtr ( „ 

nM E n daJd qtr = *(»/»(*■»■ 

This prediction is in good agreement with experiment [13]. 

3 Temperature description 

By definition, 

n 

< b ( S ) = / dcu n (q)5(qa + q b -J2v)\ A n b \ 2 , (3.1) 
J 1=1 



where is the amplitude of n creation at interaction of particles a and b. 

Considering Fourier transform of energy-momentum conservation 5-function one can 
introduce the generating function p n [14]. We may find in result that a n is defined by 
equality: 

*n(s)= r°°^e^p n (P), (3.2) 

where 

^=/{nfei)}i^- < 33 > 

The mostly probable value of (3 in is defined by equation of state: 

y/i=~\np n {0). (3.4) 

Let us consider the simplest example of noninteracting particles: 

Pn {P) = {27imK l {f3m)/f3} n , 

where K 1 is the Bessel function. Inserting this expression into (3.4) we can find that in 
the nonrelativistic case (n ~ n max ) 

a - 3 



2 (y/s — nm) 

I.e., E kin = |T, where Ekm = (y/s — nm) is the kinetic energy. 

It is important to note that the equation(3.4) have unique real rising with n and 
decreasing with s solution fl c (s,n) [15]. 

The expansion of integral (3.2) near (3 c (s,n) unavoidably gives asymptotic series with 
zero convergence radii since p n (P) is the essentially nonlinear function of (3. From physical 
point of view this means that, generally speaking, fluctuations in vicinity of f3 c (s,n) may 
be arbitrarily high and in this case f3 c (s, n) has not any physical sense. But if fluctuations 
are small (strictly speaking, they may be arbitrarily high, bu distribution in vicinity of 
/3 c (s,n) should be Gaussian), then p n {(3) should coincide with partition function of n 
particles and /3 c (s,n) may be interpreted as the inverse temperature. 

Let us define the conditions when the fluctuations are small [7]. Firstly, we should 
expand lnp„(/3 + (3 C ) over j3: 

Id 2 1 d 3 

lnp n (/3 + /3 C ) = lnp n ((3 c ) -^(3 + ^Pn(Pc) ~ ^Pn(Pc) + - (3.5) 

and, secondly, expend the exponent in the integral over, for instance, over d 3 In p n (f3 c )/df3 3 
neglecting higher decomposition terms in (3.5). In result, k-th term of the perturbation 
scries 

[ d 3 \np n (P c )/d[3 3 \ k (3k + l\ 

Pn ' k ~ i Wh^mmw I [— ) ■ (3 - 6) 



Therefore, one should assume that 

«9 3 \np n {(3 c )/d(3l « (d 2 \np n {P c )/d(3lY' 2 . (3.7) 
to neglect this term. One of possible solution of this condition is 

d 3 lnp n (P c )/d(3 3 c ^0. (3.8) 

If this condition is hold, then the fluctuations are Gaussian, but arbitrary since theirs 
value is defined by {d 2 lnp n (/3 c )/d/3 2 } 1 ^ 2 , see (3.5). 

Let us consider now (3.8) carefully. We will find computing derivatives that this 
condition means following approximate equality: 

o(3) M) (1) (M))3 

^- - r n 7 +2^-^ ^0, (3.9) 

pn pi pi 

where means the fc-th derivative. For identical particles (see definition (3.3)), 

d3g. e -/fe(9i) 



rf > (A) = „»(_i)./ |n £(9i )^__}iA 



ab\2 
n I 



/I ™ d 3 qi&~^ e ^ qi ^ I 
I n ( 27 r) 3 2g(^ ) j ■fr( 91 ' 92 '-' 9 *)' ( 3 - 10 ) 

where is the (n — k) > 0-point inclusive cross section. It coincide with /c-particle 
distribution function in the n-particle system. Therefore, l.h.s. of(3.9) is the 3-point 
correlator K 3 : 

( 3 2 3 \ 

< n^fe) >a c -3 < f[e(qi) >p c < e(q 3 ) > f3c +2f[< e( ft ) >/3c , 
i=l i=l i=l / 

(3.11) 

where the index means averaging with the Boltzmann factor exp{— (3 c e(q)}. 

In result, to have all fluctuations in vicinity of f3 c Gaussian, we should have K m ps 0, 
m > 3. But, as follows from (3.7), the set of minimal conditions looks as follows: 

K m « K 2 , m>3. (3.12) 

If experiment confirms this conditions then, independently from number of particles, the 
final state may be described by one parameter f3 c with high enough accuracy f3 c . 

Considering j3 c as physical (measurable) quantity, we are forced to assume that both 
the total energy of the system ^/s = E and conjugate to it variable f3 c may be measured 
with high accuracy 5 . 

5 Note, the uncertainty principle ~ h did not restrict AE and A/3. 



4 Real-time finite temperature generating function- 
al 



We would like to show now why and in a what conditions our ^-matrix interpretation of 
statistics is rightful. 

In modern formulations, see e.g. the textbook [16], the temperature is introduced by 
so called periodic Kubo-Martin-Schwinger (KMS) boundary condition [17]. Namely, in 
the Feynman-Kac functional integral representation of the partition function 

E{(3) = J Dfpe-Wrt. (4.1) 

the action Sp((p; z) is defined on the Matsubara imaginary time contour Cm'- (U, U — i(3), 
but fields should obey KMS boundary condition: 

<p(ti)=<p(U-i/3). (4.2) 

This is natural consequence of definition: S(/3) = Spe _/3H . 

It was offered to deform Matsubara contour in a following way: 

C m ^Csk: (ti,t f ) + (tf,ti + iP), (4.3) 

where Csk is the Mills time contour [18] and tf > U belongs to real axis [19]. Including 
the real-time parts we obtain a possibility to describe time evolution of the system 

But this attempt was not successful. First of all, we have not an evident interpretation 
ti and tf [20]. Secondly, in spite of real-time parts, this formulation unable to describe 
the time evolution [21]. 



4.1 Equilibrium media 

It was shown above that if a n is defined by (3.1) then one may introduce p n using definition 
(3.3). The Fourier transform (3.2) connects cr n p n . On other hand, p n reminds the 
partition function. 

To find complete analogy with statistical physics we should consider transition m — > n 
particles with amplitude A nm =< out; n\in; m >. Summation over n and m is assumed. 
The corresponding 5-function of energy- momentum conservation law should be written in 
the form: 

n m ,. n m 

*(Eft-Eft)=/ d4 ^( p -Eft)*( p -EPi)» p = ( E i p )- ( 4 - 4 ) 
i=i i=i J i=i i=i 

This will lead to necessity introduce independently the temperature of initial (1/A) and 
final (1//9/) states. In particle physics we can consider the final state temperature only. 

In result we get to the Fourier-Mellin transform p(/3,z) = p(0i, zf, f3f, Zf). Direct 
calculations give important factorized form: 



where the operator 

N(P,z;<j>) = J dxdx'((f> + (x)D + _(x - x',P f ,z f )(i>-(x') 



-0_(x)D_ + (x - x', ^, Zi)i+{x% i = — , (4.5) 



acts on the functional: 



p o (0 ± ) = J D ^ +D ^_ e iS^ + )-iS^-)-iV^ ++ 4> + )+iV^- + 4>-)_ (4 J3) 

At the very end of calculations one should take auxiliary variables <f>± equal to zero. 
Here D + _ D- + are the frequency correlation functions: 

D ±T (x-x',(3) = Ti J duj{q)e ±iq(x ' x ' +iT(3) z{q) 

They obey the equations: 

(d 2 + m 2 ) x G + _ = (d 2 + m 2 ) x G_ + = 0. 

So, all 'thermodynamical' information contained in the operator N(P, z; 0), but inter- 
actions are described by po (</>)■ One can say that the operator N (adiabatically) maps the 
interacting filed system on the observable states. This important property allows consider 
only 'mechanical' processes and exclude from consideration the 'thermal' ones. 

Calculating p o (0)perturbatively one can find: 

p(fi : z) = e -iV ( - ^+) +iy ( - ^-)e^-/' dxdx 'ia(x)D ab (x~x',f3,z)j b (x>) ^ (47) 
where D ++ is the Feynman (causal) Green function and 

= (D ++ y 

is the anticausal one and, as usual, j = 5/5j. At the very end one should take j = 0. 

Let us assume now that our system is a subsystem of bigger system. This would lead 
to transformation of Boltzmann factor exp{— f3e} on corresponding statistics occupation 
number N(—/3e). This means that our interacting fields system is surrounded by black 
body radiation. This is mechanical model of the thermostat (heat bath of thermodynam- 
ics). 

In result the matrix D a b takes form (we put for simplicity Zi — Zf — V): 



n(^M) n(PM)a + (Pi) 
n(P f \qo\)a-(P f ) n(^\q \ 



+2tt% -m ) ( =/o , J n _ /o N h ) (4i 



where 

a± (p) = - e !d«)l±»). 

Following Green functions: 

D ab (x - x' , 0) = J ^ e ^-^G ab (q,P) 
was introduced and the occupation number 

n ++ (qo) = n__fo,) = {e^ + ^ 2 - l}" = n(|g |^±^). (4.9) 

and 

n+-(Qo) == @ (<?o)(l + n(q p f )) + e(-g )n(-g A), (4-10) 
n_ + (g ) = e(g )n(g A) + e(-g )(l + n(-q (3 f )). (4.11) 
Assuming that $ = f3f = (3 C it is easy to find: 

G + _ (* - = G-+(t - t' - i/3) , G-+(t-t') = G+_ (t-t' + i/3) , (4. 12) 

i.e. our Green function obey KMS boundary condition. 

So, representation (4.7) with Green functions (4.8) coincide identically with (4.1) cal- 
culated perturbatively, see also [19]. 

4.2 Nonequilibrium media 

Our attempt introduce the temperature as the quantitative characteristic of whole sys- 
tem based on assumption that mean value of correlators is small. We can 'localize' this 
condition assuming that this rough description may be extended only on subdomains of 
the system. For definiteness the subdomains may be marked by space-time coordinate r. 

It should be underlined that we divide on the subdomains not the system under 
consideration but the device, where external particles are measured. Noting that external 
flow consist from noninteracting particles (including the flow of black body radiation) the 
division on subdomains can not influence on the fields interaction. 

In result we introduce the 'local' temperature l//3(r) for r-th group of interacting 
particles assuming that fluctuations in vicinity of f3(r) are Gaussian. This means that the 
mean value of correlation in the group is small, but the correlation between groups may 
be high. Nevertheless last one is not important since the external particles are on the 
mass shell. At the same time dimension of group may be arbitrary, but large then some 
r to have possibility to introduce the temperature as the collective variable. 

We can distinguish following scales. Let L q be the characteristic 4-scale of quantum 
fluctuations, L s be the scale thermodynamical fluctuations and L be the scale of subdo- 
main. It is natural to assume that L s » L » L q . 

Corresponding generating functional has the form: 



p c >i,a 2 ) = e^«p o (0 ± ). 



One may note that the 'localization' gives influence on the operator only: 

N{<Pl<p b ) = J dYdyUY + y/2)n ab (Y, y)0 6 (Y - vl% 
The occupation numbers n ab (Y, q) have same form, (3 — > (3(Y) and 

n ij (Y,y) = J du{q)e^n l3 {Y,q) 

We find calculating p perturbatively that: 

p cp {(3) = exp{-iV(-ij+)+iV(-ij-)} x 
exp{* J dYdy\j a (Y + y/2)G ab (y, {/3(Y))j b (Y - y/2)} (4.13) 

where the matrix Green function G(q, ((3(Y))) was defined in (4.8). 

5 Conclusion 

One more detail. Our consideration has show the uniqueness of Bogolyubov's solution of 
the nonequilibrium thermodynamics problem. Indeed, without vanishing of correlations 
perturbation series in the f3 c vicinity, being asymptotic, is divergent. 

We would like to stress in conclusion that Bogolyubov's creative works naturally unite 
particle and statistical physics. In result, using Bogolyubov's mathematical basis, we 
have the united scientific space in which both branches of physics, thermodynamics and 
quantum field theory, supplement each other. 
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